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Abstract. We review the apparently hydrodynamic behaviour of low transverse mo- 
mentum particles (p? < 1.5 GeV/c) produced in central and semicentral (b < 7 fm) 
heavy ion collisions at RHIC. We investigate the impact parameter dependence of var- 
ious observables, elaborating on radial and elliptic flow and particle multiplicities. We 
also discuss possible ambiguities in the initialization of the hydrodynamic system and 
Q ' present observables that should allow for their resolution. 
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c3 ' 1. Introduction and Motivation 



The first measurements at RHIC that systematically investigated the centrality dependence 
of an observable focused on elliptic flow (the anisotropic particle emission from the colli- 
sion) [[[J, followed by the centrality dependence of the absolute (charged) particle yield per 
unit of pseudorapidity [ ^ and the produced transverse energy per unit of pseudorapidity [ 
|J]. Such systematic studies of the influence of the collision centrality are of fundamental 
interest, as they represent a powerful tool to gain a detailed understanding of the collision 
dynamics: 

Firstly, non-central collisions offer additional observables due to their deformed, al- 
mond shaped overlap region, which can lead to angular dependencies (relative to the re- 
action plane) of final state observables which do not appear in central collisions with az- 
imuthal symmetry [0], Large anisotropies arise only if there is strong rescattering already 
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in the first moments ( ~ fm/c) of the collision, and (anisotropic) pressure gradients are 
building up, determining the subsequent evolution of the matter. Curiously the stronger 
forces in the direction of steepest pressure gradients lead to more transport of matter in 
those directions and thus eventually even out the differences between the radial gradients 
in the short and long direction of the initial almond. Thus anisotropies that are observable 
in the final state are built up early and in the hottest stages of the collision, as the cause of 
these anisotropies disappears during the system's evolution (on a timescale of less than ~ 4 
fm/c [ ^, H [7| ). In contrast to this self-quenching effect e.g. for elliptic flow, other dynam- 
ical quantities such as radial flow continue to grow until freeze-out and carry information 
about the full expansion stage. We explore the influence of the initial spatial anisotropies in 
terms of a hydrodynamic picture, which represents the limiting case of maximum response 
to the initially produced pressure gradients due to strong (infinite) rescattering already in 
the early stages of the expansion. Such an approach was shown to be appropriate at RHIC 
energies [ §| and is valuable to understand the global (macroscopic) characteristics of the 
expansion stage of an ultrarelativistic heavy ion collision. 

Secondly, changing the centrality leads to a varying number of participating nucleons 
and a changing size of the interaction region. The amount of energy deposited in the colli- 
sion region as well as the energy density in the system will be largest in central collisions 
and decrease with increasing impact parameter. Thus by varying the centrality, one is able 
to scan the initial energy density and in this fashion can measure excititation functions even 
without varying the beam energy. It is crucial however to disentangle such 'centrality ex- 
citation functions' from the geometric effects introduced by the varying excentricity of the 
system. In this spirit we investigate the centrality dependence of particle production per 
participating nucleon pair and transverse energy carried by the emitted hadrons to learn 
about soft and hard scattering contributions in the initial processes. 

In section |[ we introduce the underlying assumptions of hydrodynamic models and 
focus especially on different initialization scenarios. We present our results and compar- 
isons to experimental data in Section 3, which covers a discussion of particle spectra and 
radial flow, elliptic flow, multiplicities and transverse energy. Section 4 contains a brief 
summary. In the Appendix we study the effect of boost invariance on the pseudorapidity 
dependence of multiplicities and elliptic flow. This helps to understand the corresponding 
shapes of the recently presented experimental data [ ^ around midrapidity. 

2. Hydrodynamics and Initialization 

Hydrodynamics is a macroscopic approach to describe the dynamical evolution of the ex- 
pansion stage of a heavy ion collision. It is a phenomenological model that, by describing 
the evolution of thermodynamic fields like energy density, pressure, temperature and flow 
fields, circumvents the necessity of introducing unknown microscopic parameters (e.g. in- 
medium cross sections or string tensions) as required for microscopic descriptions of such 
systems. 

In the hydrodynamic description the nuclear equation of state enters the model quite 
naturally as the connection of pressure or temperature to energy and particle density. How- 
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ever it is not at all obvious that such an approach is feasible at all, as a thermodynamic treat- 
ment requires a 'large', 'macroscopic' system in local thermal equilibrium and an adiabatic 
expansion stage. But the good agreement of hydrodynamic simulations and experimental 
data from RHIC in fact point towards such rapid thermalization followed by a hydrody- 
namic expansion. 

The effective treatment however lacks a physical understanding of the underlying mi- 
croscopic processes and the early equilibration time as well as an explanation of the ob- 
viously very large rescattering rates. It must then eventually be supplemented by a mi- 
croscopic kinetic treatment to check its validity. Finally a hydrodynamic approach is only 
valid for a certain timespan of the expansion. We therefore have to introduce assumptions 
about initialization and freeze-out conditions [ |L [TT| ]. 

Existing microscopic models on the other hand so far lack rescattering mechanisms 
which are strong enough to explain the large observed anisotropies (see e.g. [ |l2[|), or the 
need to assume unrealistically large cross-sections [ Q. Only when coupled to a hydro- 
dynamic initial stage that is able to generate sufficient flow anisotropies before the system 
enters the hadronic rescattering phase, the large observed anisotropies can be recovered [ 
[l4| ]. Recently however progress in purely microscopic models was reported by introducing 
multi-Pomeron exchanges in quark gluon string models [ |5fl. 



2.1. Relativistic hydrodynamics and equation of state 

Adiabatic expansion of matter is described by the hydrodynamic equations for the conser- 
vation of energy, momentum and baryon number. In relativistic form they read 

3^ = 0, V = 0, (1) 

with the energy-momentum tensor and the baryon current 

T" v (x) = (e(x)+p(x))uf(xy(x)-gt"p(x) , f{x) = n(x)u»(x) . (2) 

These equations for the space-time evolution of the physical fields are closed by an equa- 
tion of state (EoS) relating energy and baryon density to the pressure and temperature. The 
EoS for this study contains a sharp first order phase transition, which connects a hadronic 
resonance gas at low energy densities (e < 0.45 GeV/fm 3 ) to the hard equation of state 
of an ideal ultrarelativistic gas (modeling a quark gluon plasma phase) in the high energy 
density region of the phase space diagram (e > 1 .6 GeV/fm 3 ). Further details on the equa- 
tion of state and its construction can be found in [ ^J, which also includes a discussion 
of the influence of the phase transition and details of the equation of state on final state 
observables. 

To reduce computational costs we analytically implement boost-invariance in longi- 
tudinal direction. The shape of the measured dN/dry distribution [ |l0|| indicates that 
this is well satisfied around midrapidity at RHIC energies (see Appendix). Fully three- 
dimensional calculations exist for SPS energies [ [l6[ and are under development for RHIC 
energies. 
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2.2. Initialization 



The linear scaling of particle production with the number of wounded nucleons as observed 
at SPS energies [ [l7|] indicates that the first scattering processes underlying such collisions 
are soft or non-perturbative, and that it is the number of wounded nucleons that governs 
particle production. We assume this to hold locally in the transverse plane and the number 
of produced particles to be proportional to the number of participants per unit area. For a 
collision with impact parameter b this density at a point s in the transverse plane is given 
by 



n W N(s;b) = r A (s+jl 



1-1- 



oT B (s-ib) 



B 



T B s-ib 



1-1- 



where we have introduced the nuclear thickness function 



T A (s)=J dzp A (s,z), 

with Pa(s,z) parametrizing the nuclear density profile (i.e. a Woods Saxon profile with 
appropriate parameters for a nucleus with mass number A). 

On the other hand one expects the particle production from hard collisions to dominate 
at high energies. This perturbative particle production scales with the number of collisions, 
which is given per unit area in the transverse plane by 

n B c(s;b) =cT A (s+±b) T B (s-£b) . 

In the following we study the results obtained using five different parametrizations of 
the initial state. We assume either energy or entropy density to be proportional to «wn 
(parametrizations labeled eWN and sWN, respectively) or to «bc (labeled eBC and sBC, 
respectively). In addition we use an initialization resulting from a saturation model calcu- 
lation which limits the growth of the gluonic cross-section in the transverse plane on the 
basis of geometrical arguments [jig], and label these results by 'sat'. More details on the 
initialization models can be found in [If]. 



3. Experimental Observables 

The initialization fields for energy and baryon number density are obtained by a straight- 
forward extrapolation from earlier simulations where we tuned them to fit particle spectra 
resulting from the most central Pb+Pb collisions at the highest available SPS beam energies 
[ |l9l ]. For RHIC energies, we readjust only one parameter (the initial energy density at the 
origin in b = collisions) to reproduce the final particle multiplicity observed in central 
collisions at RHIC [ ^fl]. (The equilibration time is scaled down so that its product with 
the maximum temperature is the same for RHIC and SPS systems). In our hydrodynamic 
simulations, this results in a maximum energy density eo =23 GeV/fm 3 at an equilibration 
time To = 0.6 fm/c (compared to 9 GeV/fm 3 at 0.8 fm/c for the earlier SPS simulations) 
The new initialization leads to a mean energy density of ^mean = 

4.8 GeV/fm 3 at x = 1.0 
fm/c, for which the PHENIX collaboration estimates approximately 5.0 GeV/fm 3 [ ph. 
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3.1. Spectra and radial flow 

The larger energy densities and pressure gradients in calculations for RHIC lead to stronger 
transverse expansion, which is reflected in flatter transverse mass spectra compared to SPS 
results (Fig. [l] in [ 21 1). This is in good quantitative agreement with preliminary spectra 



from the RHIC experiments [ [l(J. In the hydrodynamic simulations, the average radial 
flow velocity at freeze-out increases from 0.45 c at maximum SPS-energies to 0.55 c at 
•y/iAW = 130 GeV. The influence of the initialization on the slopes of the particle spectra 
is weak, but one observes that the spectra are getting flatter in the order eWN, sWN, eBC, 
sBC with the results for the saturation model somewhere in between sWN and eBC. This 
reflects the harder initializations of the binary collision models, which results in steeper 
initial pressure gradients and larger driving forces. 




Fig. 1. Transverse mass spectra of pions 
and antiprotons for central and semicen- 
tral collisions (the latter scaled by a factor 
0.1) as resulting from the different initial- 
izations. 



Fig. 2. Elliptic flow of charged hadrons as 
function of centrality, given by the number 
of produced particles for the different ini- 
tializations together with experimental data 




3.2. Elliptic flow 

The quantitative agreement of hydrodynamic simulations with the measured data for el- 
liptic flow, both for the momentum integrated and the minimum bias averaged differential 
elliptic flow (for not too large impact parameters and transverse momenta), was pointed out 
and discussed in an earlier work [ ||] and is reproduced in Figs . || and 9, Here we investigate 
the influence of the different initialization scenarios. 

Fig. H shows the momentum integrated elliptic flow as a function of centrality, charac- 
terized by the particle yield at midrapidity. There is good agreement with the experimental 
data for central to semicentral collisions (high to intermediate values of n c h/«max) indepen- 
dently of the underlying initialization (except of the saturated initialization, giving larger 
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Fig. 3. Elliptic flow in minimum bias 
configuration for charged hadrons together 
with the experimental data (left) and an- 
tiprotons (right). 




Fig. 4. Elliptic flow of charged hadrons as 
function of transverse momentum for fixed 
impact parameters. 



anisotropies than the others). Also Fig. |3j, which shows minimum bias averaged differential 
elliptic flow V2(pr), reveals that the results for charged hadrons are rather independent of 
the underlying initialization. Deviations can only be seen at intermediate to high pr, and 
then again especially for the saturation model. The analysis of charged hadrons is domi- 
nated by pions due to their large abundancies. Analyzing heavier particles independently, 
as e.g. done for antiprotons in the right panel of Fig. ||, shows that elliptic flow of heavier 
particles is sensitive on the details of the initialization. Fig. [5] shows differential elliptic 
flow for specific impact parameters, without avaraging over them to yield minimum bias 
results. Also here a sensitivity on the initialization is seen, especially for semicentral col- 
lisions. This sensitivity is lost when averaging over impact parameters which goes ahead 
with weighting over the resulting particles, as seen in Fig. ||. The centrality dependence of 
particle production in the different models is thus crucial and therefore studied in the next 
subsection. 

It is truly astonishing that the experimental data and the hydrodynamic results coincide 
up to impact parameters of about 7 fm and pj of about 2 GeV. Out of all models applied 
to relativistic heavy ion collisions, hydrodynamics exhibits the strongest (namely infinite) 
rescattering, and thus leads to the strongest mapping of initial coordinate anisotropies to 
final momentum space anisotropies. Hydrodynamics thus gives the upper limit of possible 
V2- That the data reach up to this limit is remarkable! 



3.3. Multiplicities and transverse energy 



Contrary to the minimum bias elliptic flow analysis, the centrality dependence of particle 
production is quite sensitive on the details of initialization. The left plot in Fig. 3.3 shows 
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the number of produced charged hadrons per participating nucleon pair as a function of 
participants from our simulations together with experimental results [ Q [loj ^Ofa - The 
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Fig. 5. Charged particle yield per participating nucleon pair (left) and transverse energy per 
produced particle (right) as functions of participating nucleons. 



experimental data disfavor the soft initialization 'eWN' and the saturation model. The 
other initializations show a tendency similar to the data, i.e. rising particle production for 
more central collisions. The data would be best described by a combination [ ^3|] of the 
two extremes for hard and soft scattering contributionsc to particle production that we have 
studied here. This is under current study within our hydrodynamic approach. 

The results for the transverse energy per produced particle are shown in the righthand 
side plot of Fig. 33. Data from the SPS [ |l7|] shows qualitative agreement with the shape 
of the curve from hydrodynamics initialized by a wounded nucleon Ansatz - a saturating 
transverse energy per particle with a wide plateau from semicentral to central collisions. 
However very recent experimental data at RHIC [ [$|] show the same shape, disfavoring the 
hard initializations of our model, which do not follow this trend. A definite statement on 
that seemingly contradictory behaviour of the centrality dependence of particle production 
and transverese energy carried per particle has to await a more careful theoretical analysis 
which is under way. 



4. Conclusions 



We have reviewed the experimental and theoretical evidence for early equilibration and 
subsequent hydrodynamic evolution of matter created in heavy ion collisions at RHIC. The 
large radial and anisotropic flow is most easily explained under the assumption of strong 
pressure gradients driving the system's expansion. The energy densities at the equilibration 
timescale reach far beyond the critical energy density. Microscopic simulations with stan- 
dard scattering cross sections fail to describe flow observables due to a lack of sufficiently 
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large rescattering. 

We further investigated a possible influence of the initialization on the observables. 
We found that no ambiguity arises for the analysis of the published results on centrality 
dependence of momentum integrated elliptic flow and transverse momentum dependence of 
elliptic flow in minimum bias configurations. On the other hand, we discussed observables 
which would allow for a distinction of the initialization, that is minimum bias elliptic flow 
for heavier particles (e.g. antiprotons), momentum dependence of elliptic flow for fixed 
intermediate impact parameters and the transverse energy per produced particle. 



Appendix 

We briefly discuss the influence of the rapidity (y) pseudorapidity (r|) transformation 
on the pseudorapidity dependence of particle yield and elliptic flow under assumption of 
boost-invariance along the beam axis. This will lead to a quantitative understanding of the 
experimental data recently presented by the PHOBOS collaboration [ p|, 

We start out with the definitions for the rapidity y := Atanh (p-/E) and the pseudora- 
pidity T| := Atanh (p z /p). From these definitions follows the connection of the differentials 
dy — ^ ■ dr\ which is the focus of interest. 

For this instructive (quantitative) excursion we make very simple model assumptions. 
We stricly assume boost invariance not only around midrapitity, but assume that all quan- 
tities are independent of y, no matter how large y is. For the spectra we use the simple 
exponentials 

dN 



.(p T ) ^Ne-V^+Pr/T 



prdpjdy 

where we use for these case studies simply T = 190 MeV and m = m n = 140 MeV. 

Boost invariant, i.e. y independent spectra transform due to the Jacobian to r| depen- 
dent spectra according to 

dN /?7-coshr| dN 

■(PT,r\) = ; = — -, — tKPt)- 



p T dp T dy\ J m 2 + p 2 cosh 2 y] PTdp T dy 



This distribution gives thus smaller values than the original one. The suppression is largest 
for low pr and vanishes for high pj. Also the reduction is larger the larger the mass is, and 
the smaller r| (i.e. around midrapidity). For large r| the distribution approach each other. 

From this it is obvious that the pj integrated spectra lead to smaller values of dN jdx\ 
around midrapidity than far away, where it approaches the constant value of dN/dy (left 
panel of Fig. ||). In reality boost invariance must break down at some rapidity, and the 
particle yield drops to zero. However the dip observed in the dN /dr\ distribution around 
midrapidity is thus well expected from a boost-invariant source, and we find in fact quanti- 
tative agreement with the preliminary STAR-data [ [l0|] from r| = —0.5 to 0.5. 

Rapidity dependent elliptic flow is defined as V2(pr,b,y) = (cos(2(|))) and the same 
for V2(pT,b,T[) where the averages are taken with respect to the dy respectively dr\ dis- 
tributions with a fixed impact parameter b. From the definitions one finds easily that 
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Fig. 6. Transformation of boost-invariant quantities (independent of rapditity y) to pseudo- 
rapidity. The left plot shows the effect of the Jacobian of the transformation on the particle 
spectra, the right plot the influence on the elliptic flow coefficient V2. 

V2(pr,b,r[) = V2(pr,b,y) = V2{pr,b,y = 0) and therefore also V2(pr,b,r[) is independent 
of the pseudorapidity. However for the total elliptic flow, when integrating over the trans- 
verse momentum, the pseudorapidity dependence will come into play, due to the different 
shapes of the spectra: 



Consider the effect of this: V2(pr) is a monotonically increasing function. For large r| 
the weighting particle distributions will become identical to the spectra dN / prdprdy and 
therefore we will recover V2 as weighted with the boost-invariant rapidity-distributions. 
Now going to smaller pseudorapidities we are weighting the differential elliptic flow with 
the suppcression of the low pj part. Thus higher p^'s with larger V2{pr) will get stressed 
and we will have larger /^-integrated elliptic flow around mid(pseudo)rapidity, than away 
from it! We illustrate this in Fig. ^, where we used for simplicity the linear relation 
v 2{ptiT\) = pr(GeV)/8.5 which is a first approximation for the experimental results on 
'minimum bias' collisions. 

Thus, a boost invariant source would show a bump in v'2(r|) and a dip in the rapidity 
distribution dN /dr\. In reality we also expect that elliptic flow drops to zero far away 
from mid-rapidity as in these regions the reaction dynamics does not allow for sufficient 
equilibration, the primary cause of elliptic flow in non-central collisions. 



V2(Tl) 
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